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1 Introduction 



The method of fundamental solutions or which is the same of discrete sources (we 
will keep to the first name) is a widely used technique for the numerical solution of 
elliptic boundary value problems which falls in the class of so called boundary methods 
reducing problems in n -dimensional domains to some equations on their (n — 1)- 
dimensional boundaries. It is applicable when a fundamental solution of the differential 
equation of the problem is known and the completeness of an infinite system of such 
fundamental solutions with singularities (sources) placed outside the domain of the 
problem is proved. The original idea of the method emerged in sixties [Q, |15|, |1(J and 
since then the method of fundamental solutions was successfully used in geophysics, 
acoustics, elasticity theory, electromagnetism and other fields. We refer the reader to 
the books [[[], |5|] and to the review |7j for bibliography and more information about the 
method. 

The aim of the present paper is to propose necessary elements for the application of 
the method of fundamental solutions to boundary value problems of electromagnetic 
scattering theory, in particular, to introduce a system of fundamental solutions for 
Maxwell's equations and to prove its completeness in appropriate functional spaces. 
Here an explanation is needed because all this seems to represent nothing new. Let 
us start with the concept of a fundamental solution for a system of partial differential 
equations with n unknowns and n equations of the form Au = , where A is a 



differential operator. Usually (see, e.g., ^FL, p. 179]) it is defined as a matrix nxn, 
denoted by $ , such that 

A$ = 5E n , (1) 

where 5 is the Dirac delta function and E n is the nxn identity matrix. Never- 
theless such a definition has no clear physical interpretation as a field generated by a 
point source, the usual meaning of the fundamental solution. In this sense the elec- 
tromagnetic fields produced by an electric and a magnetic dipoles are closer to the 
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physical meaning of a fundamental solution and sometimes they are called the funda- 
mental solutions of the Maxwell system @, Sect. 4.2], but then it is not clear how can 
they be used for the analytical solution of homogeneous and inhomogeneous Maxwell's 
equations, usually based on property ([[]). 

We propose another possibility, a fundamental solution which enjoys both prop- 
erties. It satisfies (|I|) in a sense explained below and has a clear meaning of a field 
generated by a point source. Moreover, we prove the completeness of an infinite system 
of such fundamental solutions in appropriate Sobolev spaces which makes it possible to 
apply our system to the numerical solution of boundary value problems for Maxwell's 
equations in chiral media. 

The construction of the system of fundamental solutions for the Maxwell equations 
proposed here is based on some elements of quaternionic analysis which seems to be 
the most appropriate formalism for this task. The solutions obtained are complex 
quaternions, that is instead to be a pair of three-component vectors, they have four 
components. Due to their lower singularity and simple form the numerical application 
of them is easier and more natural compared with solutions based on a matrix approach 
(see § 1). 

The main idea to obtain the quaternionic fundamental solutions for Maxwell's 
equations consists in the quaternionic diagonalization of Maxwell's equations proposed 
in [|ll|] (see also [14], [12|). The Maxwell equations for an isotropic homogeneous medium 



are reduced to a pair of quaternionic equations in which the unknown functions are 
separated. For each of these equations a fundamental solution is easily constructed and 
then linear combinations of them will give the required system of fundamental solutions 
for Maxwell's equations. The main difficulty constitutes the proof of completeness of 
this system. We base our proof on the completeness of a system of fundamental 
solutions of the Helmholtz operator in the kernel of this operator in L 2 -norm and 
make use of a quaternionic decomposition of the kernel of the Helmholtz operator. 
Our results are applied to Maxwell's equations for chiral media but they are com- 
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pletely new for a nonchiral case as well. 



2 Complex quaternions 

We shall denote by H(C) the set of complex quaternions (= biquaternions) (the letter 
EI is frequently chosen in honor of the inventor of quaternions, W. R. Hamilton). Each 



is the unit and {i k \ k = 1,2,3} are the quaternionic imaginary units, that is the 
standard basis elements possessing the following properties: 



We denote the imaginary unit in C by i as usual. By definition % commutes with 
i k , k = 0,3. 

The basic quaternionic imaginary units %\ , %i and can be identified with the 
basic coordinate vectors in a three-dimensional space. In this way a vector it from 
C 3 is identified with the complex quaternion a\i\ + a 2 i 2 + cl?,H • We will use the so 
called vector representation of complex quaternions, each a G HT(C) is represented as 
follows a = a + l? , where a® is the scalar part of a sometimes denoted as Sc(a) = ao 
and ~~ct is the vector part of a : Vec(a) = ~ct = Yl\=\ a k^k ■ Complex quaternions of 
the form a = ~ct will be called purely vectorial. 

In vector terms, the multiplication of two arbitrary complex quaternions a and b 
can be rewritten as follows: 



element a of H(C) is represented in the form a 



Yl=o a kik where {a k } C C , i 



2q — *o — ~ i-l'i i-o^k — ik^o — k — 1, 2, 3; 



1X12 = -«2«1 = ^ l 2 i3 = = H] «3«1 = -*1«3 = % 2- 



a ■ b = a b — < if, b > + if x 




+ a b + bolt, 
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where 



< if, b >:= ^^cikbk G C, 



fc=i 



if x b ] :- 



ll % 2 lz 
CL\ &2 &3 

b\ b 2 b 3 

The complex quaternion a = a — if is called the conjugate of a . Let us denote 



by (3 the set of zero divisors from H(C) . We recall (see, e.g., |Tj]) that a E & iff 
a - a = . If a ^ © U {0} then a -1 exists and a -1 = a/ (a ■ a) . 



3 Quaternionic differential operators 

We shall consider IHI(C) -valued functions depending on three variables x\ , x 2 and x% . 
On the set of such componentwise continuously different iable functions the operator 
D = Ylk=i^kdk, where dk = gf- , is defined. The expression Df , where / is an 
H(C) -valued function, can be rewritten in a vector form as follows 



Df = — div / + grad f + rot / . 

That is, Sc(Df) = — div / and Vec(Df) = grad/ + r °t / • The condition / e 
ker D is equivalent to the Moisil-Theodoresco system 

div T = 0, 

^ (2) 
grad f + rot / = 0, 

which has been studied in hundreds of works (see, e.g., 0, ||). 

Denote D a = D + al , where a is a complex constant and / is the identity 
operator. As we will see in the subsequent pages, a has the meaning of a wave 
number. Having this in mind we will assume that q^O and Im a > . 
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We have the following factorization of the Helmholtz operator 

A + a 2 = -D a D_ a = -D_ a D a , (3) 
which in particular means that any function satisfying the equation 

D a f = (4) 

or 

D_ a f = (5) 

also satisfies the Helmholtz equation 

(A + « 2 )/ = 0. (6) 
We will use the fundamental solution of the Helmholtz operator 

p ia\x\ 

0a(x) -- 



4-7T \x\ 

which fulfills the Sommerfeld radiation condition at infinity. 

Fundamental solutions K, a and /C_ Q for the operators D a and D_ a respectively 
can be obtained easily using (^). We have that the functions 

K a = —(D — a)9 a and /C_ Q = —(D + a)9 a 

satisfy the equations 

D± a JC± a = 6. 

More explicitly we have 

fC± a (x) = (±a H p — ia—)6 a (x), 

\x\ \x\ 

where x = X]fc=i x fc^fc- Note that JC± a (x) are complex quaternionic functions with 
Sc(JC± a (x)) = ±a9 a (x) and Vec(JC± a (x)) = — gr&d9 a (x) = (-^ — iaj^)9 a (x) . 
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Let us introduce the following operators 



n± Q — =F— -Dzpo 



1 

2a' 



considering them on H(C) -valued functions from ker(A + a 2 ) . Then we have the 



following statement (see the proof in |H , p. 36]). 
Proposition 1 The following relations hold: 

1. ul a = u ±a , 

2. n a n_ Q = n_ Q n a = O; 
s. n a + n_ Q = / ; 

4- as imll ±a = kerD± a we have 

ker(A + a 2 ) = ker D a © ker D- a . 

4 Quaternionic integral operators 

Let T be a closed Liapunov surface in M. 3 . The corresponding interior domain we 
denote by Q + and the exterior by Q~ . Let it be the outward with respect to Q + 
unitary normal on T in quaternionic form: it = Ylt=i n k'''k ■ Denote 

K ±a f(x) = - J K ±a {x - y)lt(y)f(y)dT y , x G M 3 \ T, 

where / is an H(C) -valued function and all the products under the integral are 



quaternionic. The following important result is well known (see, e.g., |14] , p. 70]). 
Theorem 2 Let f £ C 1 (fi+) n C(H+) and f E ker D± a (n + ) . Then 

f(x) = K ±a f(x), \/xeQ + . 
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Remark 3 In this paper the belonging of a complex quaternionic function f to some 
functional space means that each of its components fk belongs to that space. 

For the consideration of equations (f|) and (^j) in the domain Q~ one needs ap- 
propriate radiation conditions at infinity. Such conditions were introduced in Jl9| (see 



also Ul3"|). Solutions of (HI) are required to satisfy the following equality uniformly in 



all directions 



%x 1 
(1 + | — r) ■ f{x) = o(-j— r), when \x\ — > oo. (7) 



For solutions of (||) in Q the corresponding radiation condition has the form 



'1 — — — ) • f(x) = o(- — -), when Id — > oo. 
\x\ \x\ 



Then we have the following result [ 13 1 : 



Theorem 4 Let f G C 1 (r2 ) n C(Q~) , f e ker D a {Vl ) and satisfy the radiation 
condition (|7p or f e kerD_ Q ,(fi _ ) and satisfy ^). Then 

f(x) = —K a f(x) or f(x) = —K^ a f(x) respectively, for any x 6 

As we will consider H(C) -valued functions satisfying the Helmholtz equation 

(A + ct 2 )u = (9) 

in unbounded domains it will be convenient to obtain for such functions a radiation 
condition at infinity in a quaternionic form. If a solution uq of the Helmholtz equation 
is a scalar function, then the corresponding radiation condition is the well known 
Sommerfeld condition 

iau (x) — ( — — , gr&du (x)\ = o(— — ), when \x\ — > oo. (10) 

\ fI / \ x \ 

For a vector solution ~T? of the Helmholtz equation the corresponding radiation con- 
dition has the form [|, Sect. 4.2] 



. x 
rot u x - — - 

\x\ 



— — div it — ialt = °(rr when |x| — > oo. (11) 
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Let us notice that a vector solution ~Tt of the Helmholtz equation fulfills this condition 
if and only if each Cartesian component of ~Tt fulfills the Sommerfeld radiation 
condition 0, Sect. 4.2]. Thus, our quaternionic radiation condition must include both 



JIH| ) for u and ( jLT]) for Ut . 

It is easy to obtain such a condition using Proposition [I] and radiation conditions 
dj]) and (H|). From Proposition |l] we have that an H(C) -valued solution u = uq + ~vt 
of (0) has the form 

u = H a u + U- a u, 

U Q u fulfills (0) and U- a u fulfills (^). Thus we obtain 

u(x) = —- — r ■ Ii a u(x) + -^y • Ii^ a u(x) + o(- — r ), when Id — > oo. 

From the definition of H± a we have 

2au(x) = ^— ■ (D — a)u(x) + ■ (D + a)u(x) + o(-— ), when \x\ — > oo. 
\x\ \x\ \x\ 

Finally we arrive at the following radiation condition at infinity for complex quater- 
nionic solutions of the Helmholtz equation (|9|): 

x 1 
iau(x) + — ■ Du(x) = o(— — ), when \x\ — > oo. (12) 
\x\ \x\ 

As is easy to see when u = uq the scalar part of this equality gives us exactly the 
Sommerfeld condition (|10D and the vector part 



-7^- x gradu (x) 
\x\ 



o(^) (13) 
\x\ 



is a redundant equality because it is a simple consequence of the fact that a scalar 
solution of the Helmholtz equation satisfying the Sommerfeld condition at infinity can 
be represented as a single layer potential which satisfies (0). 

When u = it , the vector part of (|12"D gives us (|TT|) and the scalar part 



is again a simple consequence from the integral representation of ~Tt (see ||, Sect. 4.2] 
or 



0, P- 120]). 

Thus Qi~2"D in special cases reduces to (|10D and ( PH) and in general represents the 
radiation condition at infinity for the quaternionic Helmholtz equation. 

Note that ([/]) and (§) follow from ( |i~2|) immediately if one assumes that u G ker D a 
or u G ker -D_ a respectively. 

We will need the following operators 



S a f{x) 



-2 y K a (x-y)J?(y)f{y)dr y , xeT, 



P a = l;(I + S Q ) and Qq, = ^{1 - S a ) 

defined for example on Holder functions in the sense of the Cauchy principal value. It 
is well known that S a is a singular integral operator of the Calderon-Zygmund type 
(see || Sect. 2.5]). This implies the boundedness of the operators P a , Q a and S a 
in Sobolev spaces H S (T) for all real s . 

The following important properties of the operators P a , Q a and S a will be widely 
used in this work. 

Theorem 5 Let f G L 2 (T) . Then for an almost every point t G V there exist the 
following nontangential limits 

hm K a [f](x)=K a [f] ± (r), (14) 
and the following formulas hold 

K a [f] + (r) = P a [f}(r) = 1 -{I + S a )f(r), K a \f\-{r) = -<?«[/] (r) = - S a )f(r) 

(15) 

Let now T be a sufficiently smooth surface in order that the Sobolev space H S (T) 
for a given s be defined. 
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Remark 6 Since H S {T) C £2 CO for s > and S a is bounded in H S (T) , equalities 



JJ) and ([UJ hold for f G H S (T) , s > 



Corollary 7 1. T/ie following equalities hold on H S (T) , s > 

*5 a -^j P a Pan Qa Qai PotQot QaPa 0. 

That is P a and Q a are orthogonal projection operators on H S (T) , s > . 

2. In order for f G H S (T) , s > to be a boundary value of a function F G 
ker D a (Q + ) , the following condition is necessary and sufficient: 

f eimP a (H s (T)). (16) 

3. In order for f G H S (T) , s > to be a boundary value of a function F G 
ker D a (Q~) , satisfying (Qj at infinity, the following condition is necessary and 
sufficient: 

f eimQ a (H s (T)). (17) 

The proof of these facts in -^(r) can be found in |TJ], Chapter 5]. By Theorem [| 
and Remark |] these assertions also hold in the space H S (T) , s > . 
Needless to say that the same facts are valid for D_ a . 

5 Complete systems of fundamental solutions of 
the Helmholtz operator 

Let T be a closed surface in M 3 which is a boundary of a bounded domain Q + 
and of an unbounded domain Q~ = R 3 \ Q + . By T~ we denote a closed surface 
enclosed in Q + and enclosing the domain V and by T + a closed surface enclosing 
Q + as shown in Fig. [1]. By {Vn}^=i we denote a set of points distributed on T~ 
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Figure 1: 

and dense on T~ , and by {l/n}™=i a se t °f points distributed on T + and dense on 
T + . To each of these sets a system of fundamental solutions {0 a (x — y~)}°^ =1 or 
{0 a (x - yt)}n=i is related. Denote 6~ >n (x) = 9 a (x - y~) and 0+ n (x) = 9 a (x - y+) . 
Singularities of functions of the first system are distributed on the interior surface T~ 
and consequently every such function is a solution of the Helmholtz equation in Q~ 
satisfying the Sommerfeld radiation condition at infinity. Functions from the second 
system have their singularities on the exterior surface r + and solve the Helmholtz 
equation in Q + . 

We start with the following theorem due to V. Kupradze ]T3J, the proof can be 
found, for example, in H p. 51]. 

Theorem 8 Let T be a closed surface of class C 2 . Then the system of functions 
{9 a (x — yt)}^=i is complete in L 2 (T) . Assume additionally that a 2 is not an eigen- 
value of the Dirichlet problem in V . Then the system of functions {6 a (x — yn)}^=i 
is complete in ^(T) also. 

Our aim is to obtain a similar result for the Sobolev spaces H S (T) . This will require 
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a sequence of steps. We will show first that these systems of functions are complete 
in L 2 {VL) n ker(A + a 2 ) . Then this result will be extended to H S (VL) n ker(A + a 2 ) . 
Finally as the space H S {T) can be considered as a space of traces of corresponding 
solutions of the Helmholtz equation we will be able to prove the completeness of our 
systems of fundamental solutions for the Helmholtz operator in this space. Let us 
consider first the case of a bounded domain Vt + and then of an unbounded domain 

5.1 Interior domain 

Theorem 9 Let Q + be a bounded domain in M. 3 with a Liapunov boundary T . The 
system of functions {0a,n}^-i is complete in L 2 (fi + ) fl ker(A + a 2 ) . 

Proof. We consider the operator A + a 2 as an unbounded operator in L2(fl + ) 
with the domain H 2 (Q + ) . This operator is closed and the set L 2 (f2 + ) fl ker(A + a 2 ) 
is a subspace. Thus it is sufficient to prove that the system {0^}^^ is closed in 
L 2 (fi + ) fl ker(A + a 2 ) . Assume that there exists a non-trivial function / G L 2 (£l + ) fl 
ker(A + a 2 ) with the property 

Kn, /> i2(Q+) = for all „ G N, 

or in explicit form: 

/ 6 a {x -y+)f*(x)dx = for all n G N, 
Jn+ 

where "*" stands for the usual complex conjugation. Denote 

Vn+f(y)= 6 a {x-y)f(x)dx. 
Jn+ 

We have that Vh+/*(y+) =0 for all n G N . These equalities and the continuity of 
V n +f* imply the equality V n +f* = on T + . 

The function Vh+/* satisfies the Helmholtz equation in Q~ and fulfills the Som- 
merfeld radiation condition at infinity. Consequently Vn+f* = in Q~ . Moreover, 
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all the derivatives of Vq+/* in are equal to zero. Thus we obtain that the func- 
tion Va+f* and all its derivatives are equal to zero on T . Taking this into account 
and using the fact that V^+f* G ker(A + a* 2 ) (A + a 2 ) in Q + due to the uniqueness 



of continuation for the null solutions of this elliptic operator p0| , Theorem 6.14] we 
obtain that V^+f* = in Q + and hence / = in Q + . ■ 

Theorem 10 Under the conditions of Theorem^ the system of functions {6an}™-i 
is complete in H S (Q + ) fl ker(A + a 2 ) , s > . 

Proof. Here we use a quite general fact proved by N. Tarkhanov (for a general 
elliptic system) [22, Sect. 8.1] that a function from H S (Q + ) fl ker(A + a 2 ) belongs 
to the closure of the subspace sol(f2 + ) in H S (Q + ) consisting of all C°° solutions 
of the Helmholtz equation in a neighborhood of fi+ . That is for any function / 6 
H S (Q + ) flker(A + a 2 ) and for any e > we can find such a function /o e ker(A + a 2 ) 
in Q + , where Q + C f2 + that ||/ — /o||fl-«(n+) < s/2 . The domain Q + can be chosen 
enclosed by T + . 

For all solutions u of the Helmholtz equation in Q + we have the following estimate 
(see, e.g., [^, Theorem 11.1]) 

where the constant C does not depend on u . Due to Theorem |^ for any e\ > the 
function f can be approximated by a linear combination fx = J2n=i a n@a,n m ^ + 
in such a way that 

\\fo ~ fN\\ L2 (Q+) < e l- 

Choose E\ = e/(2C) and consider 
11/ ~ /w !!#»(«+) = 11/ - /o + /o - fN\\ H s( n+) < ||/ - fa\\ H s(n+) + Wfo ~ fN\\ H s {n+) < 



+ C\\f -f. 



N\\L 2 (a+) 



< 6 
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Theorem 11 Let T be a sufficiently smooth (the space H S {T) is defined) closed sur- 
face. The system of functions ^ s complete in H S (T) , s G R . 

Proof. For s < the result follows from Theorem |^. Let s > , given e > 
then for any u G H S (T) there exists (probably not unique) a solution of the Dirichlet 
problem (A + a 2 )U = in Q + and U | r = w , a function C/ G # s+1 / 2 (ft+) . Due 
to Theorem [10] for any £i > we can approximate it by a linear combination Un = 

^2n=l a n@ln m SUCn a Wa Y that 

||f — Un\\ hs + i/2(q+-) < E\. 

Using the continuity of the trace operator 7 we obtain that 

||w — Wjv|| HS ( r N = ||7v^ — f^Jv) Hi/s ( F ) < C \\U — C/jv||j3- s +i/2(-n+) < Ce\. 
Choosing E\ = e/C we finish the proof. ■ 

Remark 12 This theorem was proved for scalar functions from H S (T) . Nevertheless 
it is obviously valid also for H(C) -valued functions from H S (T) which in this case are 
approximated by linear combinations Yln=i c n@an where c+ are complex quaternions. 

5.2 Exterior domain 

Let Br be an arbitrary ball with a sufficiently large radius R such that Q + C Br . 
Denote = Q~ n Br . Thus is a domain in R 3 with a boundary consisting of 
T and of the sphere 8Br . 

Theorem 13 Let V be a closed Liapunov surface and a 2 be not an eigenvalue of the 
Dirichlet problem in V . The system of functions {^n}^! ^ s complete in L 2 (fi i? ) D 
ker(A + a 2 ) . 
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Proof. Assume that there exists a non-trivial function / G L 2 (tt R ) fl ker(A + a 2 ) 
with the property 



KnJ) L2{n - R) = forallnGN, 



or in explicit form: 



/ a (x - y n )f*(x)dx = for all fiGN. 



Denote 



v n-f(v)= 6 Q (x-y)f(x)dx. 

We have that V -f*(y~) = for all n e N and hence V u -f* = on T. 

The function V n -f* satisfies the Helmholtz equation in Q + . It is equal to zero 
in V due to the uniqueness of a solution of the Dirichlet problem in V and it is zero 
with all its derivatives in Q + due to the uniqueness of continuation for the solutions 
of the Helmholtz equation. Moreover in y the function V n -f* belongs to 
ker(A + a* 2 ) (A + a 2 ) . Thus due to the uniqueness of continuation for the solutions 
of this elliptic operator we obtain that V^- /* = in Q R and hence / = in . ■ 

Theorem 14 Under the conditions of the previous theorem the system of functions 
{&a,n}°°-i i s complete in H S (Q R ) fl ker(A + a 2 ) , s > . 

Proof. First, for a given e > we choose such a function f Q e ker(A + a 2 ) in 
Q R that ||/ — /q||# S (q-) < e/2 , where fl^ is a domain containing Q R and such that 

H flft = . For all solutions u of the Helmholtz equation in Q R we have the 
following estimate 

where the constant C does not depend on u . 



The proof finishes by analogy with that of Theorem 10 
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Theorem 15 Let T be a sufficiently smooth closed surface and a 2 be not an eigen- 
value of the Dirichlet problem in V. The system of functions {d^nj^-x is complete 
in H S (T) , s £l. 



The proof is completely analogous to that of Theorem 1 1 . 



6 Extensions into exterior domains 

Let us consider the exterior Dirichlet problem for the Helmholtz equation 

(A + a 2 )U = infi", 

(18) 

U\ r = u, 

and U satisfies the Sommerfeld radiation condition ( |i~0"|) at infinity. For u G H S (T) , 
s > it is known that the solution of this problem exists, is unique and belongs to a 
weighted Sobolev space in Q~ (see pT], Sect. 2.6]). For our purposes the important 



will be the fact that the solution belongs to H s+l ^ 2 (VL R ) where £l R is an intersection 
of fl~ with a ball Br of radius R chosen large enough to enclose the interior domain 
0+ . We denote by #f+ 1/2 (fi-) the union of all such H S+1 I 2 {VL R ) . 

The same will be valid if in fllTf ) we assume the functions u and U to be H(C) - 
valued and each Cartesian component of U satisfy (|TOD or which is equivalent the 
whole function U satisfy flT2"). The operator transforming u into U we denote by 



A, and as we have just seen A acts from H S {T) to H^ 1 . 

The operators Ilia introduced above act obviously from Hf jyi~) to Hj 1 ' 1 ^^) . 
Consider a function IL^Am. For u G H S (T) , s > 1 it will belong to H^ 1 ^ 2 ^^) and 
its trace (see, e.g., |2l], p. 50]) r yIl a Au G H s ~ 1 (T) . As the operators Q± a are bounded 



in H S (T) , we can introduce two new operators Q a and Q~ a as follows 
Q±a = Q±a7n±«A : H S (T) -> H s - l (T), s > 1. 
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Proposition 16 Let an H(C) -valued function u belongs to H S (T) , s > 1 . Then 

u = Q a u + Q- a u. (19) 
Proof. Consider U = Au . We have 

U = U a U + n_ a C/ (20) 

= -{K al U a U + tf_ a7 n_ a E7), 

for any point x G fi~ . Taking the limit of this equality when x tends to the boundary 
and using Theorem [5] we obtain that 

u = Qa_^Ti a A.u + Q^ a ^U_ a Au = Q a u + Q- a u. 

■ 

Remark 17 As we have seen the operators Q± a act from H S (T) to /f s_1 (r) , s > 1 , 
so equality ATQ ) can appear a little bit surprising. Nevertheless this is a reflection of the 
corresponding fact inside the domain Q~ where the differential operators H± a 

also reduce the smoothness of a function but the derivatives in ( \2dj) are cancelled. 

Proposition 18 Let f e imQ a {H s {T)) , s > . Then Q a f = Q a f . 

Proof. Let / G imQ a , that is f — Q a f . We have that the function Af satisfies 
equation (f|) in Q~ and belongs to . Moreover due to the uniqueness of 

the solution of the Dirichlet problem for the Helmholtz operator in Q~ we obtain that 
n a A/ = Af . Thus 

Af = -K al U a Af 
which on the boundary due to Theorem || gives us that / = Q a f . m 
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Let us introduce the following systems of functions 



X 



-D + a)9 a {x - y±)} 



n=l 



and 



{K\ n {x) = -{D + a)6 a {x-y±)}™ =1 



(21) 



(22) 



where the sets of points {Vn}™=i an d {y n }^ =1 are defined as in Section |[ We are 
ready to prove one of the central facts of this work. 

Theorem 19 Let a 2 be not an eigenvalue of the Dirichlet problem in V . Then the 
systems of functions {/C^ Qn }^ 1 are complete in imQ± a (H s (T)) , s > 1 respectively 
by the norm of H S ~ 1 (T) . 

Proof. Let us consider the system {fc.~ n }°° . Due to Proposition |I8| any function 
/ G imQ a (H s (T)) can be represented as follows 

/ = Qaf- 



Due to Theorem 15 for any e > there exists such a linear combination 



N 



3=1 



that 



II/-/. 



iv|lH s (r) 



< e. 



Here a~ are constant complex quaternions. Due to the boundedness of Q a we have 



/ " Qaf. 



N 



H*-i(r) 



Qaf Qaf, 



N 



_H" s -!(r) 



< C 11/ — fN\\ H s(T) < 6 



where C is a positive constant. Thus the function Q a fN approximates / in the norm 
of H'-^T). Consider 



N 



QafN = Q a f yIl 0( Ay]a j 9 a ^. 

3=1 
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It is obvious that the extension A# ■ coincides with the values of 8 a (x — jJa ) for all 
x G Q~ . We obtain 

N N N 

QJn{x) = ^2Qajtt-a(0a(x - yj))aj = — ^QalJC'^Oj = ^^2^a,j( X ) a j- 

j=l j=\ j=i 



7 Extensions into interior domains 

The results of this section and their proofs are similar to those of Section |6| and we 
present them more briefly. Here a new and natural assumption will be that a 2 is not 
an eigenvalue of the Dirichlet problem in Q + . Then for each H(C) -valued function 
u G H S {T) , s > there exists its unique Helmholtz extension, an H(C) -valued 
function U satisfying the Helmholtz equation (H) in Q + and coinciding with u on 
the boundary. As before the operator transforming u into U we denote by A . By 
analogy with the operators Q± a we introduce the operators 



P ±a = P ±al Ii ±a k : H S (T) -> H- L (T), s > 1. 

Proposition 20 Let a 2 be not an eigenvalue of the Dirichlet problem in Q + and an 
H(C) -valued function u belong to H S (T) , s > 1 . Then 

u = P a u + P^ a u. 



The proof is analogous to that of Proposition |T6 



Proposition 21 Let a 2 be not an eigenvalue of the Dirichlet problem in Q + and an 
H(C) -valued function f belong to im P a (H s (T)) , s > . Then P a f = P a f . 



The proof is analogous to that of Proposition |18[ 

Finally, by analogy with Theorem [TP] the following statement is proved 
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Theorem 22 Let a 2 be not an eigenvalue of the Dirichlet problem in Q + . Then the 
systems of functions {fC±a,n}°°=i are complete in imP± a (H s (T)) , s>l respectively 
by the norm of i/ s ~ 1 (r) . 

Remark 23 For a = a similar result can be found in jl], p. 284] ( see a ^ so references 
therein). Unfortunately the scheme of the proof proposed in that work is not applicable 
for complex quaternion valued functions due to the difficulty of introduction of an 
space which would correspond to the complex quaternionic multiplication. 

8 Complete systems for Maxwell's equations 

As we will see our approach works not only for homogeneous, isotropic, achiral media 
but also for chiral media. This last case is more general. When the chirality measure 
of a medium (3 is equal to zero we obtain the nonchiral or achiral situation. This is 
why we show our results for the case of a chiral medium, transition to a nonchiral case 
is quite easy. 

For the sake of simplicity we consider a sourceless situation. Then Maxwell's 
equations for time-harmonic electromagnetic fields in a chiral medium have the form 
(see, e.g., [0,0) 

divE(x) = div# (x) = 0, (23) 

rot E (x) = iuB (x) , (24) 

rot H (x) = -iujD (x) , (25) 
with the constitutive relations [T7J 
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D = e(E(x)+ProtE(x)), (26) 



B = H (h(x) + (3 rot H(x)Y (27) 

where u is the frequency, e and fi are complex permittivity and permeability of a 
medium and (3 is its chirality measure. 

The Maxwell equations (p3|)-(p5|) can be also written as follows 



mt E (x) = iu ^ {h(x)+ (3 mt H{x)^j, (28) 

rot H (x) = -iooe (E (x) + (3 rot E (x)) . (29) 
Introducing the notations 

E(x) = -y/ji-~3(x), (30) 
H(x) = y/e-lt(x), (31) 

we obtain the equations 

rot ijf (x) = -ia (it (x) + (3 rot it (x)\ (32) 

and 

rot It (x) = ia (~E* (x) + (3 rot ~E* (x)) , (33) 
22 



where as before a = uj^/eJI and in the case of (3 = , a is the wave number. When 
(3 is different from zero, as it will be seen below, a does not have the same physical 
meaning. There appear two wave numbers instead, a± and a-i . 

Let us consider the following purely vectorial biquaternionic functions: 



<P { x ) 



2? (x) + ilt (x) (34) 



and 



~$(x) = ~fi(x)- ilt (x) . (35) 



We have that 



Dip (x) = rot if (x) + i rot it (x) . (36) 
Using (]3~2|) and ([33]) we obtain 

Dip (x) = -(ialf (x) + a~E* (x)) - 



a(3 (d1? (x) + iDlf (x) j 



That is, 



Dip (x) = — alp (x) — af3Dlp (x) . 
Thus the complex quaternionic function ~p satisfies the following equation 

Z,+ (TT^)^ W = °- (37) 

By analogy we obtain the equation for ip 
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D - 



a 



(1 - a(5) 



if) (x) = 0. 



(3? 



Introducing the notations 



a a 



l + a/3)' ~* (l-a/3) 
we rewrite the equations ( |37l) and (|38|) in the form 



(D + ai )!?(x) = 



(39) 



and 



(D-a 2 ) if) (x) = 0. 



(40) 



When /3 = we arrive at the quaternionic form of the Maxwell equations in the 
nonchiral case, but in general the wave numbers a\ and a 2 are different and physically 
characterize the propagation of waves of opposing circular polarizations. 

As was shown in preceding sections, the functions ~Tft and if) can be approximated 
by right linear combinations of functions {^Q^.n}^ an d {^-a 2 ,n}^ = i respectively. 
The vectors ~S and ff are easily obtained from ~T$ and if) : 



^ = i(y + ?) 



and 



(41) 



2 yr T ' 2i 
Consequently, all the results of preceding sections are applicable to the electromagnetic 
field. As before we start with exterior domains. The radiation condition for the vectors 
E and H is the Silver-Miiller condition: 

, 1 



or in an equivalent form 



+ 



\x\ 



\x\ 



(42) 



\x\ 



(43) 
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uniformly for all directions. 

Note that fl42"D and ( |43l) are fulfilled automatically if as before T{? satisfies (0) and 
ip satisfies (|j). We have 

¥ = ^) = _(-_. + _. ^) + (_) 

= n ■ ^ - ^) + °(rr) = fr ■ ^ + °(rr)- 

The vector part of this equality gives us ([42]) and the scalar is a simple consequence 
of (pEHp. Starting with H instead of E we arrive at (f4"5p. 

From Theorem ||, equalities ([y]) and the last observation concerning the relation 
between the radiation conditions for ~7ft and ip from one side and for ~E* and it 
from the other, we obtain the following criterion. 

Theorem 24 Let complex vectors if and h belong to H S (T) , s > . Then in order 
for it and h to be boundary values of E* and it satisfying the Maxwell equations 
ffim and (3B[) in Q~ and fapfy at infinity, the following condition is necessary and 



sufficient 

(~e* + % h ) G imQ ai and ("e* — i h) G imQ_ a2 . (44) 

Now from Theorem [19| we obtain immediately the following important result open- 
ing the possibility to apply the systems of quaternionic fundamental solutions {^C^^j^j 
and {^-aa,™} 00 -! ^° approximation of the electromagnetic field in exterior domains. 

Theorem 25 Let both a\ and a\ be not eigenvalues of the Dirichlet problem in V . 
Then if (Kl) is fulfilled the vectors it and h belonging to H S (T) , s > 1 can be 
approximated with an arbitrary precision (in the norm of i/ s ^ 1 (r) ) by right linear 
combinations of the form 

^ N N 

3=1 3=1 
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and 

^ N N 
3=1 3=1 

where dj and bj are constant complex quaternions. 



Proof. Due to Theorem [19] we have that there exist such aj and bj that 



A? 



-? + it - 53^01^ 

3=1 



< e 



-H" s_1 (r) 



and 



N 



3=1 



< 6. 



From these two inequalities we obtain the necessary result. ■ 

In a similar way we obtain the corresponding result for interior domains. 

Theorem 26 Let both a\ and ot\ be not eigenvalues of the Dirichlet problem in 
Q + and the following condition (which is a necessary and sufficient condition of the 
extendability of the vectors if and h into Q + in such a way that their extensions 
satisfy @) and (g^|) ) be fulfilled 



-? + i h ) G mv P ai {H s {T)) and (~? - i h ) E imP_ Q2 (# s (r)), s > 1. 



Then ~e* and h can be approximated with an arbitrary precision (in the norm of 
H s ~ l (T) ) by right linear combinations of the form 

1 N N 

3=1 3=1 

and 

- N N 
3=1 3=1 

where dj and bj are constant complex quaternions. 
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9 Numerical realization 



Let T = dQ be a closed sufficiently smooth surface in M 3 . Consider the following 
exterior boundary value problem for the Maxwell equations 

rot if (a:) = -ia(lt{x) +/3rat 7?(ar)), x G f2~, (45) 
rot if (x) = ia(fi{x) +0iot~E?(x)), x G fl - , (46) 

if (a;) x rf(:r) = T(^), x G T, (47) 

where / (x) is a given tangential field and at infinity the vectors If and i? satisfy 
the Silver-Miiller radiation condition ( f4"2"D or fl4"3]). As before "rf stands for the unit 
outward normal to T. 

Using results of the preceding section this problem can be rewritten in following 
equivalent form 

(D + «i) (p(x) = 0, xeQ~, (48) 
(D-a 2 )iJ>(x) = 0, 16^, (49) 



i [{<p(x) + iP{x)) x Jt(x)] = ?(x), x G T, (50) 



Sc<p(x) = 0, Sc^(a;) = 0, a; G T, (51) 
and at infinity the functions (p(x) and ip(x) satisfy the conditions 

1 + t—/\ ■ fix) = o ( - — - ) , \x\ — > oo, (52) 



IX \ ( I 

1 — — — ) • tp{x) = O I -j — r ) , \x\ — > OO. (53) 
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Note that due to the uniqueness of the solution of the exterior Dirichlet problem for 
the Helmholtz equation condition (fTfl) implies that (p and ip will be purely vectorial 
on the whole domain Q~ . 

We look for an approximate solution of (f48|)-(|53D in the form 

N N 

<Pn(x) = ^K.-^{x)aj and ip N (x) = ^ JCZ a2ij (x)bj, (54) 

3=1 3=1 

where a 3 - and bj are constant complex quaternions. 

In order to find the coefficients a 3 - and bj we use the collocation method. In fl53|) 
we have 8N unknown complex quantities, the components of cij and bj . Thus it is 
necessary to obtain 8N linearly independent equations with respect to components 
of a,j and bj . 

Every collocation point generates four equations, two of which correspond to the 
boundary condition ( 



) and the other two correspond to (51). Consequently in or- 
der to determine the coefficients a,j and bj we need 2N collocation points. After 
having solved the corresponding system of linear algebraic equations we obtain the 
approximate solution of the problem (|45| ) - (|47|) : 



~tf(x) = —((p N (x)-if> N (x)), xeQ. 

The method described above was tested using the following exact solution. 
Let (3 = and consequently a = aci = a 2 • The vectors 



/ 



rot ~~c*6 n (x) 



c 3 d 2 9 a (x) 
cid 3 6 a (x) 
c 2 d 1 9 a (x) 



C2d 3 a (x) 



\ 



- c 3 9i6' a (x) 
-c x d 2 6 ol {x) J 
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and 

H m {x) = -— iot^ m (x), I6l 3 \{0}, 

where if G R 3 is constant, represent the electromagnetic field of a magnetic dipole 
situated at the origin j|, Sect. 4.2]. They satisfy flUf ) and ( ^6|) (for /3 = 0) as well as 
the Silver-Muller conditions at infinity. 

Let T be a unit sphere with its centre at the origin. Then and i? m give us 
the solution of the following boundary value problem 

rot i*f(a;) = —iaT?(x), x G Q~ , 
rot 7f(x) = ial^(x), x G Q~ , 



E{x) x if (x) = f (x), x G T 



where 



^ c 3 9 2 6» a (x) - c 2 d 3 9 a (x) N 



/ (x) = Cl a 3 ^(x) - c 3 <9i0 a (x) x 
y c 2 9i6> a (x) - cid 2 6 a (x) J 

As the auxiliary surface T~ containing points y~ we have chosen the sphere with 
centre at the origin and radius 0.15. In the following table we present the results 
for different values of N . The corresponding errors represent the absolute maximum 
difference between the exact and the approximate solutions at the points on the sphere 
with centre at the origin and radius 5. 
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N 


— 4 

Error for E 


— 4- 

Error for H 


3 


0.802E-02 


0.756E-02 


5 


0.346E-02 


0.250E-02 


10 


0.334E-03 


0.299E-03 


15 


0.137E-03 


0.132E-03 


20 


0.128E-04 


0.166E-04 


25 


0.187E-04 


0.149E-04 


30 


0.465E-05 


0.588E-05 


35 


0.278E-06 


0.398E-06 



A quite fast convergence of the method can be appreciated (all numerical results 
were obtained on a PC Pentium 3). 

Let us notice that the approximation by linear combinations of quaternionic fun- 
damental solutions can be applied to other classes of boundary value problems for the 
Maxwell system like for example the impedance problem. 
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